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00 . ABSTRACT 

£^ , An extensive grid of synthetic mid- and far-ultraviolet spectra for accretion disks in cata- 

clysmic variables has been presented by Wade and Hubeny (1998). In those models, the disk 
was assumed to be in steady-state, that is T e s(r) is specified completely by the mass Mwd and 
radius Rwd of the accreting white dwarf star and the mass transfer rate M which is constant 
throughout the disk. In these models, T c g(r) oc r~ 3 / 4 except as modified by a cutoff term near 
the white dwarf. 

Actual disks may vary from the steady-state prescription for T e g(r), however, e.g. owing to 
' outburst cycles in dwarf novae (M not constant with radius) or irradiation (in which case T e g in 

the outer disk is raised above T stea dy)- To show how the spectra of such disks might differ from 
the steady case, we present a study of the ultraviolet (UV) spectra of models in which power-law 
temperature profiles T c g(r) cx r~ 7 with 7 < 3/4 are specified. Otherwise, the construction of 
the models is the same as in the Wade & Hubeny grid, to allow comparison. We discuss both 
the UV spectral energy distributions and the appearance of the UV line spectra. We also briefly 



discuss the eclipse light curves of the non-standard models. Comparison of these models with 
UV observations of novalike variables suggests that better agreement may be possible with such 
modified T e g(r) profiles. 

Subject headings: accretion, accretion disks — binaries: close — novae, cataclysmic variables — stars: 
atmospheres — ultraviolet: stars 

1. Introduction the photospheric spectrum is formed. The disks in 

novalike variables (such as the prototype UX UMa 
In most types of cataclysmic variable (CV) stars and also V 603 Aql and RW Sex) should be close 

matter accretes onto the white dwarf via an accre- ^ "steady-state" where 



tion disk. In some classes of CVs, e.g. the novalike 

systems or dwarf novae in outburst, the accretion f \ -3/4 

disk can dominate the flux in the ultraviolet (UV) T e g(r) = T, 
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part of the spectrum. Several studies of the UV ^ WD 

spectra of uneclipsed CV disks have shown that 

j ,1 cjj fu where r is distance from the white dwarf and i?wn 

we do not have a confident understanding ot how . . 

is the radius of the white dwarf. The reference 
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temperature T re f depends on the mass of the white 
dwarf Mwd, its radius, and the mass accretion 
rate M: 

r 3GMwdM~ 



T re f 



87rcri? WD 



1/4 



where G and a are the usual physical constants 
(Pringle 1981). 

Modelling attempts that take into account the 
detailed vertical structure and local emitted spec- 
tra from the disk surface have so far shown con- 
siderable discrepancies between predicted and ob- 
served energy distributions or T(r) structures, es- 
pecially at short wavelength (e.g. Wade 1988; Long 
ct al. 1994; Knigge et al. 1997; Robinson, Wood, 
& Wade 1999). A general feature of the failing of 
these models is that they predict ultraviolet (UV) 
spectra that are too "blue" compared with the ob- 
servations. 

To try to improve agreement between mod- 
els and observation, the computation of the local 
emitted spectra can be made more sophisticated in 
various ways. Alternatively, additional source of 
UV radiation have been proposed to be important 
in some cases (e.g. optically thin emission from the 
base of a disk wind, Knigge et al. 1998). However, 
it is also worth considering whether the standard 
T c ff(r) formula given above should be modified. 
This may be because the disk is non-stationary 
(e.g., a dwarf nova near outburst maximum may 
not have relaxed to a steady or quasi-steady state), 
or because additional local heating besides the vis- 
cous dissipation inside the disk may alter the en- 
ergy radiated locally at the disk surface (an exam- 
ple would be strong irradiation of the outer disk 
by the central star, as considered by, e.g., Vrtilek 
ct al. 1990 in the context of low-mass X-ray bina- 
ries), or for various other reasons. 

In attempting to unravel the true structure of 
accretion disks, studies of the emergent line spec- 
trum should reveal much, because each absorp- 
tion line is formed in a relatively narrow range of 
temperatures (corresponding to a narrow range of 
radii within the disk), and thus has a character- 
istic Doppler width from projected disk rotation 
at the characteristic radius. (See Wade, Diaz & 
Hubcny 1996 for illustrative examples.) An as- 
pect of any study that considers an altered T e ff(r) 
profile should therefore be to consider its effect on 
the line spectrum from the disk. 



Wade & Hubeny (1998, hereafter WH98) pre- 
sented a grid of detailed models for accretion 
disks in CVs. The WH98 models are computed 
sclf-consistently in the plane-parallel approxima- 
tion, assuming LTE and vertical hydrostatic equi- 
librium, by solving simultaneously the radiative 
transfer, hydrostatic equilibrium, and energy bal- 
ance equations. Line transitions of elements from 
H to Ni are accounted for in the computation of 
the local spectra of individual disk annuli. Limb- 
darkening (Diaz, Wade, & Hubeny 1996), Doppler 
broadening, and blending of lines are fully ac- 
counted for in the computation of the integrated 
disk spectra. 

In this work we construct accretion disk mod- 
els for cataclysmic variables (based on the WH98 
grid and its component local spectra) , in which the 
temperature profiles T e s(r) are allowed to deviate 
in a prescribed way from the standard prescrip- 
tion 1 . We consider profiles that take a power- 
law form T e ff(r) oc r~ 7 in the outer disk, with 
7 < 3/4. This does not exhaust the possibilities 
of modified disks, but the choice is motivated by 
simplicity, and the intriguing and puzzling indica- 
tion (explored below) that such disks may in some 
sense better describe the observations of novalike 
variables. 

We discuss in the sections below our computa- 
tional method and results, including the UV spec- 
tral energy distributions, the appearance of the 
UV line spectra, and briefly the UV light curves 
of eclipsing systems. 

2. Method 

The computation of a WH98 model spectrum 
involves four stages, which we summarize briefly 
(see WH98 for the detailed discussion). First, the 
disk is divided up into concentric annuli, which 
are treated as independent plane-parallel radiat- 
ing slabs. The vertical structure of each ring is 
computed with the program TLUSDISK (Hubeny 
1990a, 1990b). Second, the program SYNSPEC 
(Hubeny, Lanz, & Jeffrey 1994) is used to solve the 
radiative transfer equation to compute the local 
rest-frame spectrum for each ring. This computa- 
tion includes typically thousands of lines selected 
from the lists of Kurucz (Kurucz 1991; Kurucz 



1 In what follows, we usually suppress the subscript on T c ff , 
but by T(r) we always mean the effective temperature. 
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Fig. 1. — The temperature profiles for a class of non-standard disks based on the the WH98 model v are 
shown. Here x is the dimensionless radius defined by x = r/RwD- The segmented appearance of the graph 
at small radii reflects the discrete set of radii used in constructing the models and is not a physical feature. 



& Bell 1995). Third, the program DISKSYN5a is 
used to combine the rest frame intensities and gen- 
erate the integrated disk spectrum. DISKSYN5a 
divides each ring into a large number of azimuthal 
sectors, and sums the individual elements with the 
appropriate area weighting and Doppler shifting 
(owing to the projected orbital motion of the gas 
within the disk). Finally, the program ROTINS 
is used to convolve the integrated disk spectrum 
with a Gaussian instrumental broadening function 
and re-sample the result uniformly in wavelength. 

In the present study, modified temperature 
profiles are generated that have T(r) oc r~ 7 in 
the outer disk, and that match smoothly onto 
the steady-state T(r) profile at a tangent point. 
The modified models therefore share the standard 
model's T(r) behavior in the inner disk and depart 
smoothly from it further out. The dimensionless 
radius of the tangent point, 

^crit = ^/-^WDi IS 



given by 

Xciit ~ ye- 8 7 ) 

(valid for 7 < 0.75). Table 1 lists x a n and other 
useful quantities for a selection of values of 7. Fig- 
ure 1 illustrates a family of modified T(r) pro- 
files for a particular model (v) from WH98. Note 
that 7 = 0.75 is just the steady-state case (as 
x cr it — > 00). We emphasize that our particular 
class of models is chosen because it allows a con- 
venient parametrization for exploratory purposes 
that has the standard model as a limiting case. 
Other classes of models can be contemplated (e.g., 
a single power-law behavior of T(r) throughout 
the entire disk) that would give qualitatively sim- 
ilar results, as well as model types that might give 
considerably different behavior. Since observed 
disks seems to differ from the standard model disks 
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by being "redder" , our considered class of models 
is possibly representative of the character of plau- 
sible changes in T(r) structure. 

We consider modifications of the standard (i.e., 
steady-state) WH98 models that have Mwd = 
1.03 Mq, namely models t, u, v, dd, and jj. The 
radius of the central white dwarf for these models 
is i?wD = 5.18 x 10 8 cm; this is also taken to be 
the inner edge of the accretion disk. For discus- 
sion of eclipse light curves in the UV, we make use 
of models based on WH98 models k and bb, which 
have M W d = 0.55 M and R WD = 9.05 x 10 8 cm. 
Table 2 summarizes some properties of these mod- 
els. For each of these standard WH98 models we 
generated modified models for 7 = 0.7, 0.6, 0.5, 
0.4, and 0.3. 

The construction of the final model spectra is 
the same as for the WH98 grid, except for the 
modified temperature profile of the disk and the 
specification of the outer radius of the disk. We 
prescribe the outer edge of the disk in two ways. 
In the first case the outer edge of the disk is 
where T eff falls below about 10,000 K, as in WH98. 
Since we consider only the far- and mid-ultraviolet 
fluxes, the restriction to temperatures greater than 
about 10,000 K results in flux errors of less than 
a few percent. In the second case, the outer 
edge of the disk is defined by the tidal trunca- 
tion radius appropriate to cataclysmic variables. 
If Mwd = 1.03 M©, the outer disk edge would be 
at about 50i?wD for an orbital period of P = 1.5 
hours or at about 100i?wD for an orbital period 
of P = 6.0 hours. (These limits were computed 
using Roche geometry and M sec = 0.1P, where P 
is measured in hours and M sec is the mass of the 
secondary star in solar units.) Table 2 shows that 
truncation of the disk at a temperature well above 
10,000 K is likely in the case of high mass transfer 
rates or low values of 7. 

Once the modified radii are computed, the in- 
tegrated disk spectrum is computed in exactly the 
same way as the standard WH98 models: local 
spectra are generated and summed together with 
appropriate area weighting, taking into account 
Dopplcr shifts (owing to orbital motions of the 
area elements within the disk) and limb darken- 
ing. We chose inclination angles of 8°, 60°, and 
76° (cosi = 0.99, 0.50, and 0.25) for illustration. 
The fluxes presented here are given for an assumed 
distance of 100 pc, and correctly account for limb 



darkening, but do not include a geometrical fore- 
shortening factor cosi. 

Note that in these exploratory calculations, we 
have deviated from requiring that the local T e g 
and photospheric gravity log g be mutually consis- 
tent: rather than compute new vertical structures 
and local spectra for each annulus, we have sim- 
ply used the vertical structures and spectra from 
WH98, shifting each model annulus to a new ra- 
dius according to its T e g. The main results of this 
study of the ultraviolet spectrum should not be 
greatly affected by this shortcut. 

3. Results 

3.1. "Color-Magnitude" Diagrams 

We show "color-magnitude" diagrams of the re- 
sulting disk spectra in Figures 2 and 3. We use 
near-monochromatic fluxes that are averages of 
f v over 10 A intervals centered at 1455 A and 
1945 A. These intervals are relatively line-free. 
Note that the flux ratio f u (U55 A)// I/ (1945 A) 
is very nearly independent of interstellar redden- 
ing (Seaton 1979). Figure 2 shows power-law 
disks that are viewed nearly face-on (cos i — 0.99) 
and are not truncated (i.e. the disks extend un- 
til an effective temperature of about 10,000 K is 
reached). The solid lines connect sequences of 
models with the same maximum temperature and 
different 7. For 7 < 0.75, the disks are hotter 
at each radius in the outer disk than in the cor- 
responding standard case, the disks extend fur- 
ther out in radius, and there is a larger propor- 
tion of the disk surface at relatively low temper- 
atures. These disks are brighter and redder, that 
is, log/„(1455A)//„(1945A) is more negative, dis- 
played further to the right in Figure 2. The dashed 
box outlines the region occupied by models when 
limb darkening is taken into account 2 . (Models 

2 To make the meaning of this region more explicit, consider 
the point marked "t" in Figure 2. This is model t in the 
case of the standard disk and for fi = 0.990. The edge of 
the dashed box descending to the right from this point is 
the track of this model as /i decreases to 0.250. The next 
vertex ascending to the right is model t with 7 = 0.3 and 
fi = 0.250. The vertex descending to the right from the 
point representing model jj (with 7 = 0.3 and /1 = 0.990) 
is model jj with 7 = 0.3 and /1 = 0.250. For all the other 
models explicitly shown (with fi = 0.990), the correspond- 
ing points for /1 = 0.500 and 0.250 lie within the dashed 
box, descending to the right from the fi = 0.990 case. 
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Fig. 2. — A "color-magnitude" diagram for the integrated disk spectra where the "color" is defined as 
log(/„(1455A)//„(1945A)) and the "magnitude" is defined as log/„(1455A) at 100 pc, with f v in mJy. The 
disks are viewed nearly face-on (cosi = 0.99) and are allowed to extend outwards until T e s ~ 10,000 K. 
The solid lines connect sequences of models with the same maximum temperature T max and different 7. The 
dashed box outlines the region occupied by models of disks inclined to the line of sight, with limb darkening 
taken into account but geometrical foreshortening ignored (see text). A spectral index a = 1/3, where a 
is defined by /„ oc v~ a , is shown by the vertical dashed line. Points representing IUE observations of the 
novalike variables RW Sex and V603 Aql are also shown (see text for details). 



with cosi = 0.50 and 0.25 are included in defin- 
ing this box. As before, the calculated fluxes are 
scaled to a distance of d = 100 pc. The purely 
geometrical foreshortening factor cos* is not in- 
cluded.) At high inclinations, limb darkening low- 
ers the disk surface brightness and reddens the 
disk. "Cooler" models such as model t and models 
with 7 = 0.3 have more pronounced limb darken- 
ing effects in the ultraviolet (see Diaz, Wade, & 
Hubcny [1996] for further discussion of limb dark- 
ening in accretion disk spectra.) 

Figure 3 is similar, with the dashed vertex at coordinates 
(—0.094, 2.717) representing the critical case for model u 
for /j, = 0.250, and the vertex at (0.105, 4.554) representing 
the truncated model jj for fj, = 0.250. 



Figure 3 shows what happens when the model 
disks are truncated at a dimcnsionless radius of 
x = 100. The dashed box again outlines the region 
when high inclination angles are considered. In 
this disk becomes brighter and redder as 7 

declines from 0.75, as long as the entire disk with 
T off (x) > 10,000 K fits within x = 100, as before. 
For 7 less than some critical value 7 cr it, however, 
the disk is no longer allowed to grow in radius, but 
each part of the disk within x — 100 continues to 
get hotter, so the disk gets brighter but bluer. At 
a fixed value of 7, truncated disks are fainter and 
bluer than non-truncated disks. 

For truncation at x = 50, corresponding to 
shorter orbital periods, the result is that the disk 
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Fig. 3. - Similar to Figure 2, except that the disks have been truncated at a dimensionless radius of x = 100. 



is smaller, fainter, and bluer than the correspond- 
ing disk truncated at x = 100. Such disks are still 
brighter than the standard disks from which they 
were derived but may be cither redder or bluer 
depending on 7 and maximum temperature T max . 

In both Figures 2 and 3, points are plotted 
which represent IUE observations of the nova- 
like variables RW Sex and V603 Aql (cf. Verbunt 
1987). The data have been corrected for extinc- 
tion and reddening and scaled to d = 100 pc and 
cosi = 0.99 (i.e., geometrical foreshortening re- 
moved) for comparison with the models. For RW 
Sex, d = 150 pc and cos-i = 0.73 are used; for 
V603 Aql, d — 376 pc and cosi = 0.96, following 
Wade (1988). The orbital periods are 5.9 hours 
and 3.3 hours, respectively. The disk models for 
the standard cases are too blue to match the ob- 
served colors of RW Sex and V603 Aql. This point 
will be discussed further below. 

WH98 discuss the contribution of the central 
white dwarf (WD) to the light from the disk (their 
Section 4.3). For the relatively hot disks that are 
considered here, the WD's contribution does not 



make a qualitative difference to the discussion, ow- 
ing to its small surface area in comparison with 
the outer disk, and has been omitted from our 
exploratory calculations. This is especially true 
for the non-standard disks of the kind considered 
here, which are more luminous than the standard 
disks. For realistic cases, the WD contribution to 
the models summarized in Figures 2 and 3 would 
be less than log f v (1455 A) w 1 or 1.5 (mJy at 100 
pc). There will be a noticeable contribution to the 
(out of eclipse) flux of a CV only when the WD is 
hotter than the inner disk, and in those cases the 
total light will be somewhat bluer than a stan- 
dard disk; in the interesting non-standard cases, 
the WD contribution is even less significant. Of 
course, the effect of the WD might be seen during 
deep eclipses of the disk, and this should be mod- 
eled in detailed studies of observed systems; this 
does not make a qualitative or important quantita- 
tive change to the conclusions of this study in any 
way. Similar considerations apply to the possible 
contribution from a boundary layer (BL) at the 
inner edge of the disk, but with even more force 
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Fig. 4. — The integrated spectra at cosi = 0.99 for the non-truncated 7 disks derived from the WH98 model 
u are shown. The fluxes are for a distance of d — 100 pc. 



because the BL radiating area would be smaller 
than the WD surface area and the UV bolometric 
correction would be very large, placing most of the 
emission in the extreme UV or soft X-ray region. 

3.2. Integrated Spectra 

Figure 4 shows integrated spectra at cosi = 
0.99 for non-truncated 7 disks derived from the 
standard WD98 model u. The disks with flat- 
ter temperature distributions (i.e. smaller 7) are 
brighter. Fluxes are again for a distance of d = 100 
pc and do not include the foreshortening factor 
cos i (which is however the same for all spectra 
shown). In addition to the overall flux level and 
color changes, specific line features that are formed 
in the outer (cooler) portions of the disk are en- 
hanced relative to features that are formed in the 
inner disk. For example, compare the C IV dou- 
blet at A w 1548, 1550 A, which is formed in the 
hotter regions of the disk, with the Si II doublet 



at A « 1527, 1533 A. In the model spectrum with 
7 = 0.7, the Si II feature is only slightly stronger 
than the C IV feature, whereas the Si II feature 
is noticeably stronger in the model spectrum with 
7 = 0.3. 

3.3. Extension to Eclipsing Systems 

Wade & Orosz (1999) discuss the modifications 
to the WH98 models to deal with the case of a 
disk that is partially eclipsed by the mass-losing 
star. We used the Wade & Orosz codes to compute 
spectra and "light curves" through various par- 
tial eclipse phases for non-standard models derived 
from the WH98 models bb (T max = 39,000 K) 
and k (T max = 29,000 K). Both models have a 
white dwarf mass of Mwd = 0.55 M Q . For the 
illustrations we choose a secondary star mass of 
M sec = 0.5 M Q , an orbital period of 0.2 days, 
an inclination of i = 71°, and an outer disk ra- 
dius of x w 51. We chose two different spectral 
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and 7 = 0.3 cases. The dashed lines in the upper plots indicate the rest wavelengths of the Si II doublet. 
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resolutions at which to sample the final model 
spectra: R = 25,000 (FWHM = 0.062 A sam- 
pled every 0.0155 A) and R = 2000 (FWHM = 
0.78 A sampled every 0.195 A). The former re- 
solving power was delivered by the G160M grat- 
ing of the Goddard High Resolution Spectrograph 
(GHRS) which flew on the Hubble Space Telescope, 
and the latter resolving power was delivered by the 
G140L GHRS grating. 

Figure 5 shows the eclipse geometry and pre- 
dicted spectra (in the restframc of the white 
dwarf) in the region of the Si II doublet (A = 
1526.71, 1533.45 A) for models bb (R = 25,000) 
and k (R = 2000) at different phases, where mid- 
eclipse is phase 0.0 (for clarity each spectrum has 
been normalized at 1600 A). We display only the 
rings in the disk that have T cS > 10, 000 K. The 
integrated disk spectra for the standard models 
(7 = 0.75) only show weak features outside of 
eclipse. As the hotter inner regions of the disk 
are eclipsed, lines that are formed in the outer 
regions of the disk (for example Si II) become 
apparent. As noted above, specific line features 
that are formed in the outer disk are enhanced 
relative to features that are formed in the inner 
disk in the models with flatter temperature pro- 
files. Hence, the Si II feature in the uneclipsed 
disks with 7 = 0.3 is relatively strong. Because 
the outer regions of the disk are relatively large in 
the 7 = 0.3 models, the change in the line spectra 
of the 7 = 0.3 models through eclipse is much 
less dramatic than in the corresponding standard 
model cases. 

Figure 6 shows continuum light curves through 
eclipse for various models based on bb and k (for 
a mass ratio Q = M\vd /Msec = 1-1 and an incli- 
nation of 71°). Here the "continuum" represents 
the average flux between 1410 A and 1530 A in 
the R = 2000 spectra. The depth of the light 
curve at phase 0.0 depends on the value of 7. The 
disk in the standard model k is totally eclipsed 
at phase 0.0 for our adopted geometry, whereas 
the flux at phase 0.0 for the 7 = 0.3 model based 
on k is roughly 35% of the out-of-eclipse level. 
Since the models with flatter temperature profiles 
have larger disks than the corresponding standard 
model, the eclipse light curves for models with low 
values of 7 are much broader than those for the 
standard models. Hence one may use the shape 
and depth of the eclipse light curve to help distin- 



guish between standard models and models with 
much flatter temperature profiles. 

4. Discussion 

There are several ways to compare the detailed 
model spectra with UV observations of disks in 
CVs. With the correct model, one should be 
able to match the overall flux, assuming the dis- 
tance and inclination of the object in question 
are well enough known. One should also be able 
to match the overall spectral shape (the "color") 
over some broad wavelength range. The predicted 
absorption line spectrum can be compared with 
high-quality observations for nearly face-on sys- 
tems or for systems in partial eclipse (Wade & 
Orosz 1999). Finally, the predicted light curves 
for eclipsing systems can be compared with ob- 
servations. In particular, the width and depth of 
the eclipse light curve are useful for distinguishing 
between various models. 

As discussed by Wade (1988), it turns out that 
UV observations of novalike variables are very dif- 
ficult or impossible to fit using steady-state disk 
models constructed from LTE stellar atmospheres. 
That is, one cannot match both the observed flux 
and color at the same time. This is still the 
case using the more self-consistent disk models of 
WH98. Indeed, this is obvious from Figures 2 and 
3 (where the open-circle symbols represent the 
steady-state disk spectra, 7 = 0.75). While the 
standard WH98 models dd and v roughly repro- 
duce the observed fluxes of RW Sex and V603 Aql, 
they fail badly in matching the colors. The present 
work suggests that "non-steady state" disks with 
a power-law T e g(r) distribution that is flatter than 
r -3/4 might fare better. For example, V603 Aql 
is on the sequence of power-law disks based on 
the WH98 model u (T max = 39,000 K) with 
0-3 < 7 < 0.4. 

There is one obvious problem associated with 
using non-standard T c g(r) profiles, either of the 
class described here or more generally, to model 
the disks in the novalike variables: what would be 
the reason that the steady-state T(r) law is not fol- 
lowed in these apparently steady-state disks? It is 
probably not reasonable to explain a non-standard 
temperature profile as being due to a mass transfer 
rate that is constant in time but varies with radius. 
This would result in sources or sinks of mass den- 
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Fig. 6. — Normalized eclipse light curves for various 7- models derived from models bb (top) and k (bottom). 



sity within the disk, in contradiction to the appar- 
ent steady-state nature of these stars. (See, how- 
ever, Knigge (1999), who discusses how a strong 
wind from the inner disk may act as a sink of mass, 
angular momentum, and energy) Hence the non- 
standard temperature profile must be ascribed to 
some other cause. One possible suspect is irradi- 
ation of the disk by a hot central white dwarf or 
boundary layer. Vrtilek et al. (1990) showed how 
irradiation can lead under certain assumptions to 
7 w 3/7 in low-mass X-ray binaries, but note that 
in that case radiation from the central source (a 
neutron star) dominates the outer disk flux, since 
the potential well into which the disk matter is 



accreted is very deep. It is far from clear that ir- 
radiation in CVs (where the central star is a white 
dwarf), acting alone, will be strong enough to en- 
force such a flat temperature profile, and we do not 
claim to have shown that irradiation is the cause of 
a flattened T(r) profile 3 . However, if irradiation is 
the cause of the non-standard temperature profile, 

Indeed, we do not even make the claim, on the basis of 
the limited work described here, that a flattened T(r) pro- 
file is the basis of the discrepancy between standard disk 
models and observations of luminous CVs. We have only 
demonstrated consistency of such modified disks with ob- 
servations of two test objects, in terms of UV colors and 
fluxes, along with pointing to other more sensitive tests 
involving spectra and eclipse light curves. 
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then our models are incomplete — we have merely 
assigned the same local spectra (computed under 
the assumption of no irradiation) to different radii 
in the disk. A correct treatment of irradiation 
would include its effect in the vertical structure 
of the disk (cf. Hubeny 1990a) and the emerg- 
ing spectrum would differ from merely a sum of 
hotter "normal" spectra. An irradiated disk spec- 
trum might even resemble to some extent a sum 
of Planck spectra, which Wade (1988) showed did 
match colors and fluxes of novalike variables in 
the ultraviolet spectral range. Until detailed spec- 
trum models incorporating a correct description 
of irradiation are developed, the "problem of the 
novalike variables" must be considered to still be 
unsolved. 

5. Summary 

We have presented a class of model spectra of 
accretion disks which have non-standard power- 
law temperature profiles (i.e., T e fj(r) oc r~ 7 with 
7 < 0.75). When 7 < 0.75, the disks are brighter 
and redder than disks in the corresponding stan- 
dard case, assuming the disks with the flatter tem- 
perature profiles extend outward until a tempera- 
ture near 10,000 K is reached. Non-standard disks 
that are truncated as some smaller radius are still 
brighter than the standard disks from which they 
were derived but may be either redder or bluer 
depending on 7 and the maximum temperature 
T max . For smaller values of 7 (flatter temperature 
profiles), specific line features that are formed in 
the outer disk are enhanced in the integrated spec- 
trum relative to features that are formed in the 
inner disk. For eclipsing systems we have illus- 
trated how the eclipse light curves for these non- 
standard disks are broader and less deep than the 
light curves for the corresponding standard disks. 
Such patterns may be useful in diagnosing whether 
disks in CVs do in fact have non-standard T(r) 
laws, or whether the failure of standard model 
disks to fit the data lies elsewhere. 

Finally, we point out that detailed standard 
disk models (T e fj(r) cx r~ 3 / 4 ) of Wade & Hubeny 
(1998) cannot simultaneously match the observed 
UV fluxes and colors of the novalike variables 
RW Sex and V603 Aql. Wade (1988) reached a 
similar conclusion, using disk spectra synthesized 
from stellar atmosphere models. Our present ex- 



ploratory work suggests that a disk with a much 
flatter temperature profile, T e s(r) oc r~ 7 with 
7 « 0.35, provides a better match to both the flux 
and color of V603 Aql. However, the WH98 mod- 
els as applied in this case are incomplete; since no 
mechanism is identified that results in the flatter 
temperature profile, we have simply constructed 
models as if there were a mass transfer rate that 
varies with radius (an apparently unphysical situ- 
ation for the novalike variables). If the flat tem- 
perature profiles which seem to be favored by this 
exploratory analysis are in fact due to irradiation 
(although for energetic reasons this seems prob- 
lematic), then further refinements to the models 
should include a detailed and consistent treatment 
of irradiation. 
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Table 1 

Data for Power-Law T(r) Disk Models 



7 


^crit 




T(x = 50)/T max 


T(x = 100)/T max 


0.7 


12.250 


0.288 


0.108 


0.066 


0.6 


3.361 


0.678 


0.134 


0.088 


0.5 


2.250 


0.848 


0.180 


0.127 


0.4 


1.842 


0.929 


0.248 


0.188 


0.3 


1.633 


0.969 


0.347 


0.282 



Note. — See text for definition of a; cr i t . The dimensionless radius is 
x = r/Ry/D, where r is the radial coordinate in the disk. T max is the 
maximum effective temperature of a steady-state disk, attained at x = 
1.36. 



Table 2 
Properties of Disk Models 



Model 




^rnax 




T(x = 50) a 


T(x = 100) a 


Family 


logM 


(K) 


x(T = 10,000iT) 


(K) 


(K) 


t 


-10.0 


29,330 


9.40 


10,180 


(8270) 


u 


-9.5 


39,110 


14.30 


13,580 


11,030 


V 


-9.0 


52,160 


21.70 


18,100 


14,710 


dd 


-8.5 


69,560 


28.80 


24,140 


19,620 


jj 


-8.0 


92,760 


43.70 


32,200 


26,160 


k 


-9.0 


29,330 


9.40 


10,180 


(8270) 


bb 


-8.5 


39,110 


14.30 


13,580 


11,030 



Note. — For model sequences t, u, v, dd, and jj, Mwd = 1-03M Q and 
]?wd = 5.18 x 10 s cm. For model sequences k and bb, Mwd = 0.55M Q and 
-Rwd = 9.05 x 10 8 cm. The steady-state mass transfer rate M is in units of 
solar masses per year, x = r/i?wD, where r is the radial coordinate in the 
disk. T max = T c h(x = 1.36). All temperatures are effective temperatures. 
See WH88 for further details. 

a T{x = 50) and T(x = 100) are shown for the case 7 = 0.3. Models 
discussed in the text are extended to T c s = 10, 000 K or arc truncated at 
x = 50 or x = 100, if T c ff = 10, 000 K has not been reached at that radius. 
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